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Abstract 

Using soft-collinear effective theory (SCET), we examine the 1/rrif, corrections to 
the factorization formulas for inclusive semi-leptonic B decays in the endpoint 
region, where the hadronic final state consists of a single jet. At tree level, we 
find a new contribution from four-quark operators that was previously assumed 
absent. Beyond tree level many sub-leading shape- functions are needed to correctly 
describe the decay process. 
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1 Introduction 



Inclusive semi-leptonic and radiative B meson decays offer many opportunities to test 
the flavour sector of the Standard Model. The radiative decay B — > can be used as 
a probe for new physics, and the most precise measurements of the CKM matrix element 
\V u b\ are based on the semi-leptonic decay B —> X u iv. However, experimental cuts on the 
final state are needed to suppress a strong background from charm production. These 
cuts often force the kinematics of the decay into the so-called shape-function region, 
where the hadronic final state is collimated into a single jet, which carries a large energy 
of order m&, but a small invariant mass of order mbAQCD- While the conventional operator 
product expansion (OPE) breaks down in this region of phase-space, it is still possible 
to make rigorous predictions using a twist expansion, which sums singular terms in the 
OPE into non-local operators evaluated on the light-cone PQI21EI- This method amounts 
to using QCD factorization formulas which separate the physics from the disparate mass 
scales m\ 3> i^Aqcd Aq CD into hard, jet-, and shape-functions, respectively. To the 
above mentioned decays these ideas have been applied in [3J E] . The recent development 
of the soft-collinear effective theory (SCET) as applied to inclusive decays jHJ 13 El M 
has provided a natural framework from which to prove these QCD factorization formulas 
to all orders in perturbation theory, and also to sum the Sudakov logarithms appearing 
therein; such studies have been performed to leading order in l/m& |7| ITU} ITT]. 

The goal of the present work is to use SCET to analyze the structure of the factoriza- 
tion formulas beyond leading order in l/m^. The motivation for this is two-fold. First, 
from the phenomenological side, the inclusion of power-suppressed effects will soon be 
necessary to keep theoretical predictions on par with the improving experimental accu- 
racy being achieved at the B factories. Second, from a purely theoretical standpoint, 
although the machinery needed to begin a discussion of power corrections to QCD fac- 
torization formulas within SCET has been available for some time, there have been no 
theoretical efforts which apply this to a concrete example. For inclusive decays a study 
of sub-leading shape-function effects was carried out previously in O HHJ H3J [J3] by 
more elementary methods, and restricted to the tree approximation. In this paper, we 
take the first steps towards analyzing power-suppressed effects in inclusive B decays in 
the shape-function region beyond this approximation by an investigation of the general 
structure of factorization, and an enumeration of all relevant shape-functions. We per- 
form a complete calculation of the hadronic tensor only at tree level. However, even 
in this approximation we find a new effect related to certain four-quark operators, and 
some discrepancies with earlier results. 

The factorization formulas are derived by performing a two-step matching from 
QCD — ► SCET — ► HQET, identifying the hard functions in the first step of matching, 
and the jet- and shape-functions in the second. The first matching is done at scales of 
order m&, the second at scales yW>A~Qci> Since the hard and jet-functions can therefore 
be computed perturbatively in the strong coupling, it is evident that factorization ex- 
tends to sub-leading order in 1/m.f,. The advantage of the effective theory framework is 
that it provides a transparent book-keeping of all the relevant interactions at every step 
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in the calculation, including power-suppressed effects, and allows us to identify all the 
HQET operators that define the sub-leading shape-functions with no restriction to tree 
level in the hard (/i ~ m&) and hard-collinear (// ~ ^?77^Aqcd) fluctuations. 

The content of the paper is as follows: Sections EJ to 0] provide an outline of the 
structure of factorization at sub-leading order in l/m bl including a basis for the SCET 
currents and HQET shape-functions needed at order l/m b) but to any order in the 
strong coupling a s . In Sectional we pick one of the many terms in the expansion and 
demonstrate its reduction to a convolution of a jet- and a shape-function. In SectionElwe 
compute the structure functions that parameterize the semi-leptonic and radiative decay 
spectra in the tree approximation at order l/m b . We find many structural simplifications 
at tree level, but also a tetra-local term that has been missed before. Sections [7| and |H1 
follow up on these observations by commenting on the nature of possible simplifications 
beyond tree level and by providing a numerical estimate of the correction from the new 
tetra-local contribution. We conclude in Section El 



2 Factorization at sub-leading order 

The effects of QCD in semi-leptonic B decays are contained in the hadronic tensor W^ v . 
We can calculate the hadronic tensor via the optical theorem by taking the imaginary 
part of the forward scattering amplitude, which we define as 

W" v = -Im (B(v)\T^\ B(v)). (1) 

7T 

We use the state normalization (B(v)\B(v))=l and drop the velocity label from now on. 
The correlator T^ u is the time-ordered product of two flavour- changing weak currents 

= ij d 4 xe- iq - x T{J^{x), J"(0)}, (2) 

where q is the momentum carried by the outgoing leptonic pair (B — > X u lv) or photon 
[B — ► X s j). The calculation of this correlator is the essential element to obtaining the 
decay amplitudes. Any differential decay distribution can be derived from the compo- 
nents of the hadronic tensor. We shall use the notation of ^H] and write the hadronic 
tensor in terms of five scalar structure functions 

W/u, = Wi {p^Vy + v^v - 9^vp - ie^ af3 p a v p ) (3) 
- W 2 g^v + W 3 v ll v v + W 4 {ppv,, + v^pv) + W 5 p^p u , 

where the independent vectors are chosen to be v , the velocity of the B meson, and 
p = rrifyV — q with m& the b quark pole mass. (We use the convention e 0123 = — 1.) The 
Wi are regarded as functions of p 2 and vp. The relation of these partonic variables to 
the final state hadronic invariant mass P 2 and energy vP is 

vp = vP - (M B -m b ), p 2 = P 2 - 2 (M B - m b )vP + (M B - m b ) 2 . (4) 
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Of course, there is a dependence of the Wi on the current J assumed in (J2J), which is 
not made explicit in our notation. Since we are interested in either the semi-leptonic or 
radiative decay, the relevant currents are 1 



The decay kinematics are assumed to be in the so-called shape-function or SCET 
region, where the hadronic jet emitted at the weak vertex carries a large energy of 
order rrib, but a small invariant mass squared of order ra&AQCD- Such a jet is referred 
to as hard-collinear, and implies the existence of three widely separated mass scales 
m 2 ^> tb&Aqcd Aq CD , along with a small parameter A 2 = Aqcd/^6- The soft-collinear 
effective theory (SCET) offers a tool with which to calculate the decay amplitudes as a 
series in this small parameter A. An arbitrary momentum p is decomposed as 



where n± are light-like vectors satisfying n + n_ = 2. For inclusive decays it is sufficient 
to use what is referred to in the literature as SCETj, which contains only hard-collinear 
and soft degrees of freedom. The components of a hard-collinear momentum are defined 
to scale as (n + ph c , Phc±, n-Phc) ~ A, A 2 ) and those of a soft momentum as p s ~ 

m(,(A 2 , A 2 , A 2 ). From now on we will refer to hard-collinear momenta as simply collinear, 
and to SCETi as SCET. We find it convenient to work in a frame of reference where 
v j_ = and n_t? = n + v = 1. Furthermore, the SCET expansion in A refers to a frame in 
which the total transverse momentum of the hadronic final state is at most of order A. 
For any given q, we therefore choose the frame where q± = [rribV — p}± = —p± = and 
compute the invariant components of the hadronic tensor in this frame. 

It has been shown that the decay amplitude factorizes into a convolution of hard, 
jet, and soft functions at leading order in A jHHj. Using SCET, our aim is to show that 
an analogous factorization holds at sub-leading order in A. The scaling properties of the 
collinear jet imply that no Lorentz invariant quantity can be formed at order A, so the 
leading power corrections appear at order A 2 . In technical terms, transverse momenta of 
order A appear only in the internal integrations over the collinear momenta in a jet, and 
since the integral over an odd number of transverse momenta either vanishes or must be 
proportional to one of the external soft transverse momenta of order A 2 , the resultant 
expansion is in powers of A 2 ~ 1/ttj&. We therefore need to go to second order in the 
SCET expansion. Before giving any explicit formulae, we will outline a procedure that 
as a matter of principle could be used to establish factorization at any order in A. 

1. Matching to SCET/HQET. In the first step we remove the hard scale TTtu clS db 
dynamical scale by matching the QCD Lagrangian and currents onto their corresponding 

For radiative decays the restriction to currents implies that we consider only the contribution from 
the leading electromagnetic penguin operator in the effective weak Hamiltonian. The complete result 
contains additional terms from four-quark operators. 







(6) 
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expressions in HQET and SCET, where fluctuations are characterized by the jet scale 
^Aqcd- We then use these effective theory quantities to calculate the correlator in (J2J. 
In the following we use the position space formulation of SCET [HUH], which is especially 
well suited for the study of power corrections, since the power-suppressed Lagrangians 
and currents are already known. 

We first discuss the SCET/HQET Lagrangian. Its explicit form to order A 2 is [9 j 

C = £ [in-D 4 ^ C -^W) ^~\tr (Fr^) 
+ h v ivD s h v + q s iip s q s - i tr (FfF^) 

-|- -|- -|- L YM T 2^ + 1 ^YM ' ^HQET- I ' ) 

2=1 

The £ denotes the collinear quark field, q s the soft quark field, and h v the heavy quark 
field of HQET. The covariant derivatives are defined as iD c = id + gA c and analogously 
for iD s , but a -D without subscript contains both the collinear and soft gluon field. The 
quantities F* {F°) are the field-strength tensors built from the soft (collinear) gauge 
fields in the usual way, except for the definition of F°, where n^D rather than n-D c 
appears. The collinear and soft fields are evaluated at x, but in products of soft and 
collinear fields the soft fields are evaluated at xt = (n+x/2) n_ = x + n^L, according to 
the multipole expansion. In this notation scalar, while xt is a vector. 

The power-suppressed terms in the effective Lagrangian read 

Zxt^WcgF^W^, 
l^xn^_W c gF^W^, 

q s WliipU-li9± c W c q s , (8) 

(2) 

where the W c are collinear Wilson lines. We have omitted the terms CY q , since their 
field content implies that they do not contribute to the current correlator at order A 2 . 
The explicit expressions of the Yang-Mills Lagrangians Cy^ can be found in They 
are needed only in the calculation of l/m^ corrections beyond tree level. The SCET 



C 



(2) 

H 



c 



(1) 



4 



Lagrangian is exact to all orders in perturbation theory, receiving no radiative corrections 
jHj. In contrast, the A 2 HQET Lagrangian is 



C 



(2) 

HQET 



2171b 



h v (iD s ) 2 h v + %^ Ka^gFfhv 



(9) 



where C mag 7^ 1 represents the renormalization of the chromomagnetic interaction by 
hard quantum fluctuations. 

A second source of hard corrections is related to the matching onto the SCET heavy- 
light currents. The QCD currents Jj = ipTiQ are represented in SCET as convolutions 
of dimensionless short-distance Wilson coefficients depending on quantities at the hard 
scale m 2 with current operators J,- composed of HQET and SCET fields. The matrix 
elements of these effective currents are characterized by fluctuations on the order of the 
jet scale m fe AQ C D and the soft scale Aq CD . We write this convolution as 

0/T 4 Q)(x) = e- im " v - x £ C§\Si, ■ ■ ■ , B n ) ® jf \h, ■ ■ ■ , s„; x), (10) 

where the (g) stands for a convolution over a set of dimensionless variables §i = SiVfib- The 
superscript k refers to the scaling of the current operator with A relative to the leading- 
power currents, and the subscript j enumerates the effective currents at a given order in 
A. As with the SCET Lagrangian, the collinear fields in the current operator J- depend 
on x + Sin + but the soft fields including h v are multipole-expanded and depend only on 
x_. The factorization formula will eventually be formulated in terms of convolutions 
over longitudinal collinear momentum fractions, which we define as Ui = n+pi/rrib. The 
momentum space coefficient functions are related to those defined above by 

C^\u l )= [Yld8 i C$\8 i )e5 HUt . (11) 

i 

A basis for the most general set of order A 2 currents including radiative corrections 
has not yet been discussed in the literature. We shall return to this point in Section 
|3] For now, we simply note that these short- distance Wilson coefficients depend only 
on quantities at the hard scale m 2 and are identified with the hard functions in the 
factorization formula. 

We have now achieved the factorization of the hard scale from the soft and collinear 
degrees of freedom and can write the correlator as 

T>» = H jf (h, ■ ■ ■ , §n) ® Tp^i-h, . . . , s n ), (12) 



The hard function Hjji is a product of SCET Wilson coefficients C\j'C\,y (along with a 



'ij ^i'j' 

possible C mag from the HQET Lagrangian) and T^, ,fJ-l/ is a correlator of SCET currents. 
In what follows we will suppress all indices and denote the correlator by T eff . 
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2. Collinear-soft factorization. In the second step we factorize the matrix element of 
T eff into soft and collinear pieces. Towards this end, we first redefine the collinear fields 
according to [7j 

£ = Y^°\ A C = YA®Y\ W C = YW^Y\ (13) 

and immediately drop the superscript on the redefined fields. The Y above is a soft 
Wilson line involving n_A s evaluated at X-, and the redefinition of the collinear Wilson 
line follows from that of the gluon field. The effect of this redefinition on the SCET 
Lagrangian and currents is to transform every n_D into n_D c , and to replace q s (h v ) by 
Y^q s (Y^h v ) and iD s ^ (F s ^) by YHD^Y (Y^F° u Y). Note that the collinear Wilson line 
appears only in and W}iD c ^W c , so the positions of soft and collinear Wilson lines 
in any field product can always be inferred from the transformation of the fields under 
the collinear and soft gauge symmetries. 

We shall treat the power-suppressed Lagrangian terms in the interaction picture. 
Since after the collinear field redefinition the leading order SCET Lagrangian (see (J7J)) 
no longer couples the collinear to the soft fields, we can factorize the matrix element. 
The B meson by definition contains no collinear degrees of freedom, meaning that the 
B meson state is represented as the tensor product \B) ® |0), where the first (second) 
factor refers to the soft (collinear) Hilbert space. It follows that the matrix element of any 
SCET current correlator, including in general time-ordered products with sub-leading 
interactions from the Lagrangian, can be written as 

(B\T cS (s 1 ,...,Sn)\B) =i J d 4 xd*y . . .e iimiV -^ x (B\h v [so&Mds]h v \B)(x^, ?/_,...) 

x (0 1 [collinear fields] 1 0)(si, . . . ,s n ;x,y, ...). (14) 

The additional integrals over d 4 y . . . are related to insertions of the power suppressed 
Lagrangian. The soft matrix element depends only on n + x,n + y, . . ., so the integrations 
over transverse positions and the n_ components can be lumped into the definition of 
the collinear factor. The soft and collinear matrix elements are then linked by multiple 
convolutions over the light-cone variables x + ,y + , . . .. 

3. Definition of "shape-functions" and u j et-f unctions" . The matrix element of the soft 
fields between B meson states defines non-perturbative (leading and sub-leading) shape- 
functions S(n + x, . . .). We define these in momentum space according to 

S(x 1+ ,...,x n+ ) = J duj 1 ...du n e-^ + - M ^S(u 1 ,...,uj n ). (15) 

We shall see below that at order A 2 we can have up to a triple convolution over the 
variables u>i, coming from two insertions of the order A SCET Lagrangian. These scalar 
functions need to be modeled and introduce hadronic uncertainties into any phenomeno- 
logical applications. Properly identifying the independent functions is therefore an im- 
portant task to which we devote Sections 0] and [7| 
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The matrix elements of the collinear fields between the collinear vacuum involve 
fluctuations with virtualities mbAq CD . After integration over position arguments they 
define perturbatively calculable jet-functions, which have the form 

3(u 1 ,...,u i ;u 1 ,...u n ). (16) 

(Alternatively, we shall use the variables p 2 i2 ., where p Ul2 . = p — u 12 ...jn + /2 and 
W12...3 = oj\ + uJi + ...ujj.) Evaluating these functions removes the collinear degrees 
of freedom and defines the final step of matching, SCET — > HQET. The complicated 
form of the power-suppressed SCET currents and Lagrangian imply that many new jet- 
functions appear at order A 2 and we will not list the general set in this paper. Instead, we 
simply note that they are perturbative objects and so introduce no inherent theoretical 
uncertainties. 

Having carried out these steps, we can express the correlator in terms of a factor- 
ization formula. We insert the definitions of the hard coefficients, jet-functions, and 
shape-functions and work in momentum space. A generic term in the factorization for- 
mula is a sum of convolutions over hard, jet, and shape-functions 

T = ^2H(ui, ...Ui)® 3(Ui, . . . ,Ui]UJi, . . . ,UJ n ) ® S(ut, • • • ,^n)- (17) 

To discuss more precisely the structure of this factorization formula, we need to identify 
the set of jet- and shape-functions appearing at order A 2 . For the jet-functions, the 
necessary step is to derive the full set of heavy-light currents in the presence of radiative 
corrections. We turn to this topic in the next section. 

4. Caveats to factorization. The operators whose matrix elements define the jet- and 
shape-functions have singularities related to the light-cone expansion. Above it has been 
assumed that there exists a regularization and subtraction procedure that is compati- 
ble with the properties of the Lagrangian crucial to factorization. To the best of our 
knowledge, dimensional regularization is adequate for this purpose, but we do not have 
a general proof of this statement. 

The factorization formula (JT7j) is composed of convolutions over the soft light-cone 
variables cjj, and we must assume that the convolutions of the perturbative jet-functions 
with the shape-functions converge. Little is known about the functional dependence of 
sub-leading shape-functions, but a divergence of the convolution for Ui — > would be 
surprising, since it would indicate that long-distance physics is not accounted for by the 
HQET Lagrangian. On the other hand, divergences for u>i — > oo, if they existed, would 
be of short-distance nature and could presumably be treated with a modification of the 
factorization procedure. 

3 SCET currents 

In this section we discuss the matching of the QCD heavy-light currents onto SCET. The 
order A operators have been investigated in several places [HI 13 \U\ including one-loop 
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radiative corrections to their coefficients [TBI EE ED], but the order A 2 case has been 
given only at tree level [HI E] ■ To give a general factorization formula for the correlator, 
we will need a complete basis for the heavy-light currents in the presence of radiative 
corrections. 

The QCD currents Jj = ibTiQ are represented in SCET as a convolution of di- 

(k) 

mensionless short-distance Wilson coefficients with a current operator J- composed of 
heavy-quark and SCET fields, as in (fTUj). Our task is to find a set of current operators 
jj up to order A 2 relative to the leading-order currents. We work in the frame where 
v± = 0, which considerably reduces the number of allowed currents at a given order in 
A. 

In currents at position x, all soft fields are evaluated at x_ = (n + x/2)n_, but the 
collinear fields may be shifted along the light-cone to positions x + Sin + owing to the 
non-locality of SCET. The variables s» are integrated in a convolution product of coeffi- 
cient functions and current operators, which allows us to eliminate factors of in + d and 
l/(in + d) from the operator by a redefinition of the coefficient function. In the following 
we use the convention that operators do not contain in + d (or its inverse) operating on 
collinear fields, and denote the position argument x + Sin + of collinear field products by a 
subscript "s" . It is also convenient to write down the operator basis independent of the 
Dirac structure Tj of the QCD current. In this case we assume that the short-distance 
coefficients are Lorentz-tensors CV,-..., where the dots stand for further indices on the ef- 
fective current. We can always decompose these into scalar functions using nt, i/ 4 , g^ v ', 
and i& vpu . The scalar functions depend of course on the Dirac structure of the QCD 
current. 

Leading order. At leading order the only possible SCET currents are 



Here Tj = {1, 75, 7 a± } denotes a basis of the four independent Dirac matrices between £ 
and h v . For a given QCD current, we would proceed to decompose the tensor coefficient 
functions into scalar ones. For instance, for the V-A current relevant to semi-leptonic 
decay the three tensor coefficient functions C M , C 5fl , and C mi ^ multiplying the jf^ may be 
expressed in terms of three scalar coefficient functions multiplying the current operators 



(18) 




(19) 



Next-to-leading order. At relative order A the basis consists of the currents 



(i) 



(£W c ) s r j (x ± D a h v ) 



1 



4 1} = {IW^std^TjK 

4 1] = {lW c ) Sl [WliD^ c W c ] S2 T 3 h, 



(20) 
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The above currents should be multiplied by an overall factor l/m& to make the coefficient 
functions dimensionless. The basis of operators with scalar coefficient functions can be 
found in [TT1 H%1 IT9"] . Here we use a basis where the transverse derivative is taken outside 
the collinear Wilson line. This has the advantage that only the third operator has a tree 
level matrix element with a transverse collinear gluon. 

The first two operators are "two-body" , that is, they depend only on x_ and one other 
position x + sn + . The coefficient functions of any two-body operator to any order in the 
A expansion can be related to the coefficient function of the leading-power currents. The 
reason for this is that they all follow from matching the QCD matrix element (q\ Jj|6) to 
SCET, and the number of independent form factors needed to parameterize the QCD 
matrix element is equal to the number of independent currents at leading power |18j . 
The coefficients of the two-body operators follow from the expansion of the full QCD 
form factors. The operator Jg is a three-body operator, whose coefficient is determined 
by the expansion of (qg\Ji\b) with a transverse gluon to leading power. 

Next-to-next-to-leading order. Many structures are possible at this order. We find 



Ji 



(2 



J. 



(2 



J. 



(2 



7(2 
J i 

7(2 



(2 



J- 



(2 



7(2 
J S 

7(2 
J 9 



J 



7(2 
J ll 

7(2 
J 12 

7(2 
J 13 



7(2 
J 14 



2 (£W c ) s r i (n_xn+.DA) 
^(£W e ) s r j (x ±li x lj ,DZD v 8 h v ) 

(tw c ).id'' ± r j (x ± D a h v ) 

(CWcjsid^id^Tjhv 
{£W e ) ai \y^iDLW c ] aa T j {x ± D a h 

(£W c ) ai [id»(W^Dl c W c ) S2 
(£W C ) S1 [[W^iD^WcU [W!iDt c W ( 



CJS3 



r h 
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(21) 



The trace refers to colour, and an overall factor 1/ra 2 should be included to make the 
coefficient functions dimensionless. The argument for the completeness of this basis is 
as follows. The QCD matrix element (q\Ji\b) is expanded to order A 2 , which involves 
the transverse momentum squared, n_ times the collinear momentum, or the heavy 
quark residual momentum. Together with the terms from the multipole expansion of 
the heavy quark field, this accounts for the two-body operators j{ 2 \. Next, we can have 
(qg\Ji\b) expanded to first order in the transverse momentum of the collinear quark or 
the collinear gluon, which gives Jj 2 \, with coming from the multipole expansion. 
The coefficient of Jj 2 ^ is related to that of the order A three-body operator . Some 
of the scalar coefficients of Jg 2 g are also related to However, contrary to the case 
of two-body operators it is no longer true that the coefficients of all sub-leading three- 
body operators are related to the leading one, because the number of invariant form 
factors needed to parameterize (qg\Ji\b) in QCD is larger than the number of the leading 
three-body currents. We also get new three-body operators from (qg\Ji\b) with a gluon 
corresponding to a n^A c field or n-A s , but gauge invariance requires that the soft gluon 

(2) (2) 

is part of a covariant derivative, so we get the three-body operator J 10 . In fact, J5 

(2) 

and J 10 could be eliminated using the collinear quark and gluon equations of motion. 
Finally, we get four-body operators from (qgg\ Ji\b) with two transverse collinear gluons 

or from (qqq\Ji\b) with three collinear ( ant i-) quarks. There are several possible colour 

(2) 

structures, which can be chosen as in Jn_ 15 . 

It is useful to gain some intuition as to how this general set of heavy-light currents 
translates into terms in the factorization formula. First, it is obvious that insertions of 
heavy-light currents which contain no soft gluon fields simply build up a set of power- 
suppressed jet-functions convoluted with the leading order shape-function. This includes 
two insertions of or single insertions of J^8-i5- ^ we wor k at tree level, only two- 
body operators can contribute to the current correlator, and we need only consider two 
insertions of J^l or a single insertion of j[ 2 \. The full set of three- and four-body 
operators is only needed when one aims at an accuracy a s AQ CD /mft in the calculation of 
the hadronic tensor. 

4 Basis of shape-functions 

In this section we will collect the results for the independent matrix elements of soft 
fields ( "shape- functions" ) needed to parameterize the hadronic tensor at order A 2 . The 
possible time-ordered products that build up the current correlator to this order are 



a 



+ sym. 
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c) jWjWcW +sym. 

e) JWJWCWCW. (22) 

Inspection of the effective currents and Lagrangian shows that from these products we 
obtain the following soft operators (leaving out colour and spinor indices): 

1. From a) and b) 

(h v Y)(Y^h v ), (h v Y)(YUD»h v ), {Ki-ipy^K), 
{h v {-ip" s \-ip V ;Y){Y^h v ). (23) 

We leave out operators related by hermitian conjugation. The convention is such 
that fields in different parentheses stand at different positions 0,x-, Z-, . . . and 
that colour indices in brackets are contracted. Hence all these terms are bi-local. 

2. From the single Lagrangian insertions c) and d) in addition 

(h v Y)(Ykgn»_F° ± „Y)(Ylh v ), {h v {-i)^Y){Y^gn»_F^ v Y){Y^h v ), 
(h v Y) (YUgn^F^Y) (Y%), (h v Y) (YUgF* ±l/± Y) (Y*h v ), 
(h v Y)(Yl[iD?\ign»_F° ±u }Y)(Ylh v ), (h v Y)(Y^in.D s ,igF^ ± ]Y)(Y^h v ), 

{h v Y){Y^h v )iJd'zC^{z). (24) 

These are tri-local (with exception of the last line). The second-to-last operator 
comes only from the Yang-Mills Lagrangian £ym- 

3. From the double insertions e) 

(h v Y)(YHgn^ ±v Y)(YHgn^ ±a Y)(Y^h v ), 
(h v Y)(Y^ v )(q s Y)(Y^q s ). (25) 
These are tetra-local. 

We note immediately the increase in complexity of the matrix elements needed to pa- 
rameterize the 1/mb corrections to all orders in perturbation theory, since there appear 
tetra-local light-cone operators including a four-quark operator which has not yet been 
discussed in the literature. 

It is convenient to parameterize the shape-functions in a covariant derivative basis. 
Reinstating colour and spinor indices, the field products listed above can be derived from 

(h v Y)(x_) aa (Y^h v )(0) bl3 , 

Ch v Y)(x^ aa (YUD^Y)(z^ cd (Y^h v )(0) bf3 , 

(^y)(x_) aQ (y t i^ ± y)(z 1 _) cd (y t ^ ± y)(z 2 _) e/ (y t ^)(o)^, (26) 
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mil 



and the expressions in the last lines of 
identities Y'in-D.Y = in_d, 



YUgn u _F° u Y = - [m_<9, YUDfX 



respectively, using in particular the 



in 



-dY^[iD s Y) 



(27) 



The derivatives are meant to act on everything to the right irrespective of position 
argument, but derivatives in square brackets (not to be confused with commutators) as 
in the last expression act only within the square brackets. With this convention 



(YHD%Y)(z) cd = id»5 cd + (Y^{tD»Y})(z) cd . 



(28) 



Note that since in the second term the derivative acts only on Y, this term is a colour 
octet (in light-cone gauge, n_A s = 0, and we have Yr[i.D£Y] = gA£). The first term is a 
colour-singlet, but does not depend on z. Hence if, for instance, (YHD^Y)(z) cd appears 
in a tri-local matrix element, the colour- singlet part is only bi-local. 

We decompose the matrix elements of (|2*Bj) into a set of scalar shape-functions. The 
only possible Dirac structures between two static quark fields h v . . .h v are 1 and (7^ — 
y^)75- Defining e^ u = ie^ pa n p _v cr , we write (recall that xt = x + nt) 



(B\Ch v Y)(x_) 



N c 2 



S(x+), 



/3a 



(29) 



(B\(h v Y)(x_) aa (Yih v )(0) bf3 i I d 4 zC ( mET (z)\B) 



1 il/H 



2m b N c 2\ 2 



HQETV 

Skin( a; +) + Cmag(^fe/A t )5 ma g(3 ; +)] 



'Pa 



(30) 



(B\CKY)(x-) aa (YUD^Y)(z^ cd (Y^h v )(0)^\B) 

1 / 1 + j \ \ s ba s cd 



2 \ 2 L 1 N c 

\ / Pa K L 



+ 



O rpArpA 

z 1 ba 1 cd 



N? - 1 



iS'(x+)v^+ (iS'(x+) -T 1 (x + ,0))n 



T 1 (x + , z+) n - 



+ 2 \ 2 



-7px7s- 



2 R 2 

I Pa 



+ 



O rpArpA 

L 1 ba 1 cd ~ 



T 2 (x + ,z+) }, 



(31) 



(s|(^y)(x_) aa (yti^ ± y)( Zl _)^(yt i ^y)( Z2 _) e/ (yt /iu )( )^| j e) 

1 A + V X 



2 V 2 2 

\ / Pa 



" I (•'' • ) 
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+ 



+ 



+ 



N*-l 2 

2 SefT^T^ 1 p ~ 

iV c 2 - 1 2 
N? — 1 



C/l(x+, 2 2+ ) - C/ 2 (^+, ^2+) - 2 2+ , 2 2+ ) 

C/i(x+,2;i + ) + U 2 (x+,z 1+ ) -U 1 (x+,z 1+ ,zi+) 

I4l(x+, Zi+, Z2+) — U2(x + , Zi+, Z2+) — Us(x + , Zi+, Z2+) 



Aif ABC T A T B T C 4Ncd ABC T A T B T C 

■U 2 {X + , Z 1+ , Z 2+ ) + ~—z — W 3 (£ + , Z 1+ , Z 2+ J 



iV c (iV 2 - 1) 



+ - 



^1,2,3 — * ^4,5,6 f , 



/3a 



(iV?-l)(JV?-4) 
i | Ml o, tf 1>2 - C/ 3 ,4, 



(32) 



(s|(^y)(x_) aQ (yt^)(o) 6/3 (5 s y)(z 1 _) C7 (yt ?s )( Z2 _) d(5 | J B) = 

) 13a 



iV c 2 [2 V 2 



5a 7 Vi(x+, 2 2+ ) + y«5 7 V 2 (a; + , zi+, 2 2+ ) 
+ i- 5 -y ^3(x+,z l+ ,z 2+ ) + (r/-y)5 7 V 4 (a;+,^i+,22+) 

(7^75)57 Vs( x +, z i+, z 2+ ) 
+ (^7^5)57 V 6 (x + , z 1+ , z 2+ ) + (^-7^75)57 V 7 (x+, z 1+ , z 2+ ) 
+ Vs(x+, zi+, z 2+ )] \ 



1 (i + i u i + y 



/3a 



+ 



4 T^TQ. 1 l + f 
(N*-l)\2\ 2 , 



'/3a 



V1-4 



9-12 



/3a 



V, 



13-16 



(33) 



This decomposition makes no assumption on whether the external state is a B meson, 
a B* meson or a 6-hadron, but it assumes that it is averaged over polarizations, so 
that the only available vectors are v and n_. Our notation uses lowercase letters for 
shape-functions depending on a single variable, capital letters for those depending on 
two variables, and calligraphic letters for those depending on three variables. The only 
exception is our notation for the leading order shape-function S(x + ), where we use a 
capital letter even though it depends only on a single variable. The multi-locality of a 
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given shape-function is determined by decomposing the colour structure of the covariant 
derivatives (YHD%Y)(z)cd as in (J28|) . and then using that the singlet component has a 
lower degree of non- locality than the octet component. The parameterization is chosen 
such that the colour contractions with the tree-level jet-functions, as well as the limits 
z = in and z\ = Z2 in ()32j) . take a simple form. For instance, writing the 
combination t(x+) — T 2 (x + , 0) for the colour singlet term in the fourth line of (jHljl ensures 
that only t(x + ) appears at tree-level. We used the heavy quark equation of motion to 
reduce the number of independent functions. The gluon field equation can be used to 
eliminate (Y*in + D B Y)(z) a i- This does not lead to a simplification in practice, and we 
have therefore kept T\( X-\- , ) clS db basis shape-function. 5" is the derivative of S with 
respect to x + . 

It follows that to order 1/mb, but to arbitrary order in a s in the coefficient functions, 
the differential decay rates depend on a large number of multi-local shape-functions. 
Fortunately, the structure of the result is much simpler in the tree approximation as we 
shall see below. In addition to displaying the general set of shape-functions, the above 
enumeration of soft matrix elements allows us to clarify in more detail the structure of 
the convolutions in the factorization formula ()17j) . Since the effective currents contain 
at most three products of collinear fields at different positions, the number of collinear 
convolutions is at most two. The shape-functions are at most tetra-local, so the maximal 
number of soft convolution integrals at order l/m& is three. This results in the structure 

T = H ■ 3(u) <8> S{u) 

+ H(ui, u 2 ) <8> 3(ui, u 2 ; uj) <g> S{uS) + ^ H( u ) ® 3(u; uu 1} u 2 ) <g> S(ui, u 2 ) 

+ #(^1,^2,^3) ® S(ui,u 2 ,U3,) + • • • , (34) 

where the ellipses denote 1/mf terms not considered here, and for each term the most 
complicated structure is shown. The momentum space shape-functions are defined in 
terms of the coordinate space shape-functions given above by the Fourier transform (|15|) . 
The variable u)i corresponds to n_ki, where ki is the (outgoing) momentum of soft fields. 

5 Example of factorization at order l/ra& 

We now have all the ingredients needed to calculate the factorization formula to order 
A 2 ~ l/m&. To illustrate the step- by-step procedure outlined in Section |2J we work out 
the contribution from the time-ordered product 

r{jW f (x)J 2 (1) (0)z / d 4 zCf\z)} (35) 

as an example. All other terms can be treated in a similar way. However, there is a large 
number of them, and we do not list them in this paper explicitly. 
Specifically, we consider the currents 

</ (0) = {£W c ) si rh v , J 2 (1) = {ZW c ) sl i*$ ± rK, (36) 
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which appear in the SCET expansion of the vector current. In the first step we obtain 
the corresponding contribution to the current correlation function T pu (p = m^v — q) , 



T" v = i J d 4 xe vpx J dhd^C^is^C^ih) 

T{h v {xJ)^{Wli){x + s x n + ) [(£W c )(s 2 n+)iQ d 



l u h v (0) 



(37) 



The coefficient function can be related to as explained in Section El but the 
detailed form of this relation is not important for the following. We now use translation 
invariance to shift all fields by the amount —s 2 n + , then perform the change of variables 
x — ► x + (s2 — s\)n + , z — > z + S2n + . This does not affect the position of the soft fields, 
since [sjn + ]_ = 0. The integrations over §1^ can then be performed to obtain 

= H(n + p/m b )i J d 4 x e ip ' x T{h v {xJ)^{Wtt){x) [(£W c )(0)if ± ]Y K(0) 

i J d'z m c )( Z )zi^_gF^z^(W^)} (38) 

with H(n + p/m b ) = (n+p/ m^C^ ( n +P/ m b) ■ {C^ is real.) The effects of the hard 
scale have turned into a multiplicative factor H rather than a convolution, because both 
effective currents were only two-body. Next we perform the collinear field redefinition 
that decouples collinear and soft fields. We also use (|27)l and integrate by parts to let 
the m_9( 2 ) act on the collinear fields. The result is 

(B\T^\B) = H(n + p/m b )i J d A xd*z e ip ' x z pl - 

(B\(h v Y)(x^) aa (Y^[iD s p Y])(z-) cd (Y^h v )(0)^\B) 



(39) 



We now insert ()31|) for the B meson matrix element of the soft fields. Only the colour 
component T^T/J contributes, and we are left with the term parameterized by the tri- 
local shape-function T 2 . Going to momentum space shape-functions, we find 



(B\T^\B) 



N2-1 2 2 



" I 7 -^^7 7s7 M 



H(n + p/m h ) 



/3a 



J d^d^n^i,^) J d*xd 4 Z e iP-* e -^+ e -i^+(-i) z P± 



(40) 



in. 



A')(v\t{[(ZW c )C9j]{o^ 
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The integral over x and z defines the momentum space jet function fl. In the frame with 
p_L = it can depend on the external momenta n + p, p 2 , and the convolution variables 
£Ji )2 , and the integral equals 

^Y^{^ Y± ) ftn+prfw,^). (41) 
Inserting this into (jjQ)) we obtain the final result 

H{n + p/m b ) J duJxduJ 2 T 2 {uj 1: uj 2 ) d(n+p,p 2 ;ui, u 2 ), (42) 

which expresses the contribution from (|H5|) to the current correlator as a product of a 
hard coefficient and a two-fold convolution of a shape-function and a jet-function. 

It is instructive to examine the result in the tree approximation for the jet-function. 
In this approximation 

ntvee, 2 \ lU+pin+p U +P U +P 

d {n+p,p ;wi,w 2 ) =w 2 — 2 5— = -5 5~ > ( 43 ) 

where we used n + puo 2 = p 2 ^ - pl 12 (p Wl ... n = p - + . . . + cu n )n + /2). Each of the two 
terms after the second equality depends only on a single combination of the ooi j2 , which 
allows us to write 

J duJiduj 2 T 2 (uj l: uj 2 ) d trce (n + p,p 2 ; uj u uj 2 ) 

= J du^f J duj'(T 2 (LU-uj',uj')-T 2 (uj,uj')). (44) 

The structure of this result implies in coordinate space that rather than the tri-local 
matrix element (J3T|) parameterized by a function of two variables T 2 (x + , z + ), the tree 
level approximation requires only the simpler combination T 2 (x + , x + ) — T 2 (x + , 0) related 
to the bi-local matrix element 

(B\Ch v (-i)lllY)(x„) lu MY f h v )mB) - (B\{h v Y)(x^) lu MY ] ^h v m\B). (45) 

The simplification of the tree level result is due to an equation of motion identity behind 
the cancellation of propagators in (jlHJ)- We can see this directly in coordinate space by 
using the SCET equation of motion for the free-field propagator. Defining the contraction 

Z(x) aa Z(y) bl3 = iA(x-y)5 ab ( 1 ^-) , (46) 

II \ A / af3 
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Figure 1: Tree diagrams contributing to the current correlator T^ u . Not shown are 
diagrams that vanish when n + A c = 0, n_A s = 0, or are symmetric to those shown. 

the function A (z) satisfies 

in-dA(z) = 5 {4) (z) - ^¥-A(z). (47) 

in+o 

When this is used on the product of the two collinear propagators in the second 
term on the right hand side of ()47|) gives zero, because we work in the frame p± = 0, and 
the two delta-functions produce the two terms T 2 (x + , x + ) — T 2 (x + , 0). We will see more 
of these simplifications when we work out the complete result in the tree approximation 
in the following section. In Section [7| we ask in more generality whether the appearance 
of tetra- and tri-local shape-functions in the formalism could be spurious, and whether, 
perhaps, the final result could be expressed in terms of only bi- or tri-local matrix 
elements beyond the tree approximation. 

6 Tree approximation 

In this section we calculate the current correlator T^ u (J2J) and the hadronic tensor in the 
tree approximation including the l/m^ power corrections. In this approximation we can 
set the collinear gluon field A c to zero except in £>fg, where we can draw a tree graph 
with external soft light quarks (see Figure [IJ. The weak current Jj = ipTiQ including 
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the hard coefficient functions (at tree level) is given up to order A 2 by [HI El 

Ji = j {0) + E + E 4 2) (48) 



i=l i=l 



with 



1 ^ 



d 2 



l lb Vl 



J? } = CT," n + D s h v , J 2 (2) = fcZ^DL D» ±s h, 

4° = -Zi^T^^T^D^hv, J 4 (2) = (49) 

We now compute all the time-ordered products including insertions of sub-leading 
Lagrangians that follow from the generic expressions (j22j) and the explicit form of the 
Lagrangian (JSHSJ) and the currents The diagrammatic representation of this com- 

putation is shown in Figure The final result is rather compact, but the intermediate 
expressions among which many cancellations take place are lengthy. To display this 
result, we introduce a short-hand notation 

J ] a J b iC c = i J d 4 xe ipx T^jl(x)J b (0)i J d i z£ c (z)} (50) 

and similarly for the other terms. We also define the integral operators (x + = (n + x)/2, 
P± — 0, n+p > 0) 

/ 2 -/(, + ,.-/Ae-/|le-^-/( I+ , 

= i dx + e m '^f(x + ), (51) 
Jo 

1J , 



i ^ 



1*2— I dz + f(x + ,z + ), (52) 
n + p Jo 



2 The last diagram in Figure^is not reproduced by a correlation function of effective currents in SCET 
but requires an additional local term that we have not discussed. Since this term has no discontinuity 
it does not contribute to the hadronic tensor and decay distributions. 
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and present the intermediate result in light-cone gauge n_A s = (Y s = 1). With this 
notation, the leading-power result for T^ u reads in the tree approximation 

j(o)t j(o) = j* h v ( x _)rl^Tjh v (0). (53) 

Following the classification of time-ordered products in (j22j) . we calculate the indi- 
vidual contributions and obtain 

(I) = jf )f J(°) = 0, 

(ii) = jf = -/* 1 (^(-^ s± ) 2 )(x_)rl^r^„(o), 

(in) = jft j(o) = _/* JL ^hSI 1 ))^)^^^, 

TI-l.jP Z Z 

(iv) = jfVw + jwtjf = JLi* ^(-i^M^l^Mo) 

(V) = jPjW = -/* — (^(-^: ± ))(x_)ri^ 7 ^r,/ i ,(o), 

(vi) = j^jP = 0, 

(vii) = j^j^hcf = -/*2(^(-^: x ))(x_)rlt rj ^ ± (^)^(o) 
+ /* A foH^)^' j(x_)r}^rA,(o), 

(vni) = (jp + jwt jW) = /* ^((^^)(x_)ri^ 7 ^r,^(o) 

-^(^-)ri^7^r,(^)(o)), 
(ix) = j(°)tj(°)i4V = / 2 *r{/ax_)rttrA(o)i / c^4 2 <Ww}, 

(X) = J(°)tj(°)j£g) = 0, 

(xi) = jWjWirtg = -/^^^^(x^rJ^r.-^^n^F^j^)^^), 

(xii) = j(°)t z4f = /* I ^(x_)ri^ 7 ^ 7 , ± ^F;^ ± (,_)r J ^(o), 

(XIII) = i J(°)t J(°) zX^ z£f) = - J(°)t J(°) icg> 
z 
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+ i* 2 (h v9 A^)(x„)rl^-r jg A; ± (z.)h v (o) 

-I* 1 {KgA^gAlMx^i-TM®) 



k 2 



-i;h v (x^T]^T,(gA^gA^ - [id^gA; ± ])(z-)h v (0) 

+ 1* 1 (^(^^ - [^^j))(a;_)r.i^r,^(o), 



(XIV) = J(°>t J^iC^iC^ 



d 4 k d 4 h d 4 k 2 
(2vr) 4 (2tt) 4 (2tt) 4 



d 4 xd 4 Zl d 4 z 2 e ipx J 



(54) 



n + p jo 



In order to arrive at this result we have made repeated use of the equation of motion 
f!47jl for the collinear propagator. The expressions simplify considerably after adding 
(III) + (V) + (VIII) and (II) + (VII) + (XI) + (XII) + (XIII). Together with the other non- 
vanishing terms (IV), (IX) and (XIV) we obtain for the current correlation function 
at order l/m&, 



il/mb 



- tree 



dx+d-U i(n_p-Lo)x + \_ 



2tt 



uj + ie 



x 



r{(^y)(x_)ri^r i (yt /lu) (o) % J d 4 z c$ QET ( 



+ 



2m^ 
1 

n + p 



(^(-i^jy)(x_)rttr i (yt^ )(0 ) + (^Y)(x_)rt^r,.(rt^A)(o) 
;(^(-i^ ± )Y)(x_)rt^ 7 , x ^r i (yt /lv)(0) 

+ (h v Y)(x_)T\^ 1 ^T J (Y\^ ± )h v )(0) 



4- 



+ 



n + p 



+ 



£ + ^ + (^y)(x_)rtt(yt ( _^ s±)( _^ g jy )(z _ )r . ( yt^ )(0) 

- [ X+ dz 1+ f X+ dz 2+ g 2 (h v Y)(x^)Tt r ^Y ± T A (Yk)(z 2 -) 
+P Jo Jz 1+ 2 



n + p 



x (gV)(^ 1 _)T^ t± ^r j (Vt^)(o) 



(55) 
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We have reinserted the soft Wilson lines Y, which makes this expression gauge invari- 
ant. The derivatives are understood to act on everything to their right (or to their left 
when indicated by the arrow) independent of the position argument of the field. It is 
worth noting that the double insertion of Cjp (XIII) was expected to be tetra-local. 
However, only the delta-function terms survive in the application of the equation of 
motion identity and the complexity of this term is reduced to a tri-local term of the 
form h v (x-)[iD±iD±](z-)h v (0). On the other hand, the four-quark contribution from 
the double insertion of jQ* is tetra-local. 

We have also performed the calculation in a general frame where p± 7^ 0, in which the 
result is a significantly more complicated expression. In particular, the double insertion 
of gives a tetra-local term proportional to p\. Superficially, this seems to require 
a much larger set of shape-functions, including a different degree of locality. However, 
the hadronic tensor depends only on two kinematic invariants vp and p 2 , so two out 
of the three variables n + p, n_p and p\ must be sufficient to reconstruct the complete 
information. We therefore conclude that the specific convolutions of jet-functions and 
tetra-local shape-functions from the double insertion of (and other similar terms not 
present for p± = 0) in the general frame cannot contain independent non-perturbative 
information despite their appearance. We can see this technically by noting that the 
transverse momentum is defined with respect to a choice of vectors n_, v . A transverse 
Lorentz boost from a frame with p± = to a frame with p± ~ A can be effected by a 
reparameterization of n_ — > n_ + 2e± — e\n + with e± ~ —p±/n + p and v fixed. The 
complete SCET expansion is invariant, because it reproduces a Lorentz-invariant theory, 
but the transformation reshuffles terms in the A expansion. In particular, the leading- 
order Lagrangian changes — > + e±5C^ + . . .. Accordingly for matrix elements 

(O) -> (O) + e ± J d 4 z (T{Oi6£(°\z)}) + ..., (56) 

from which it is clear that the leading-power bi-local term in the frame p± = gives 
rise to tetra-local terms proportional to p\ at order A 2 in the general frame. While this 
explains the structure of terms in the general frame, we did not verify explicitly the 
equivalence of the expressions, which seems to be a technical but unilluminating task. 

We now proceed to the evaluation of the hadronic tensor (JIJ. Starting from (J55|) this 
requires that we i) take the imaginary part, ii) take the B meson matrix element and 
iii) insert the decomposition of the soft matrix elements into scalar shape-functions as 
given in Section 0] The imaginary part of is obtained by the replacement 

^— - -7T%). (57) 

uj + te 

The B meson matrix elements can all be expressed in terms of the previously defined 
shape-functions, if the relation 

(B\(h v Y) (x-) aa {YHD^Y) (zS) cd (Y ] h v ) (0)bp\B) 

= (S|(^y)(-x_)^(F t (-z)^F)(^ - x^) dc (Y ] h v )(<d) aa \BY (58) 
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and d fl± (B \ . . . \B) = is used. The scalar decomposition contains odd terms propor- 
tional to e^, which can be eliminated using 6^7^75 = ~ We then 
express the position space shape-functions in terms of their Fourier transforms defined 
as in (JT5j) and perform the integrations over the positions left in . The result for the 
hadronic tensor including now the leading-power contribution (in the tree approxima- 
tion) becomes 



= J du 5{n_p - u) I 



-tr 
2 



1 



P+Vl—Tj iS(w) + [s kin (^) + C ma g(m fe //i)s mag (^)] 



u s (m) + 2v s (uj) 
n + p 



1 

+ tr 

4m& 



- ^rjtrJ - t(u)) + — — - — t ]■ 



+ 



2n + p 



-tr 



4 



p +7 ^7 5 rj 7a± 7 5 r, t(u) 



P+^-75-P+rJ— 75^ (« fl (w) - 2u«(o;)) 



(59) 



where P + = (1 + ^)/2, and we introduced the definitions 

u s (lu) = / duoiduo 2 J 2 (co; u} 1} to 2 ) {ux{oJi)8{u 2 ) + U 1 (cu 1 , co 2 )), 



u a (u) = J duidu 2 J2(u;wi,V2)U 3 (uJi,u} 2 ), 

v s (u) = / dujiduj 2 du z J^ux,^,^) g 2 V 10 (uJi,LJ 2 ,LJ 3 ), 



v a (uj) = / duJiduj 2 du 3 J 3 (uj-,uj 1 ,u 3 ,uj 2 )g 2 Vi3(Lj 1 ,uj 2 ,Lj 3 ), 



(60) 



J 2 {n-p\uj u uj 2 ) = -—1m ^ +P } = — (<J(n_p-wi - w 2 ) - <J(ra_p - wi)), 

w x 1 r (^+P) 3 

J 3 (n_p; wi,o;2,a;3j = — Im ^— 2" — 2 — 

5(n-p — LUx) 5{n_p — uji — co> 2 ) 5{n_p — uj\ — uj 2 — uj 3 ) 
u 2 (uj 2 + uj 3 ) uj 2 uj 3 ^3(^2 + ^3) 

-n 2 5(n^p - ^1)5(^2)^(^3), (61) 

and used that t(a>) is real ^2j- The denominators in the definition of J 2 and J3 are 
understood to be supplied with a principal value definition. Eq. (f5*H|) is our final result, 
valid for arbitrary Dirac structures Tj and Tj of the weak currents. Despite the fact that 
the hadronic matrix elements are tri- and tetra-local the final result can be written as 
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single convolutions of integrated shape-functions. Hence, from a phenomenological point 
of view, the l/m& corrections can be parameterized by a set of functions depending only 
on a single variable, just as at leading power. We discuss in Section [7| whether this holds 
when loop corrections are included. The factor g 2 in the definition of the integrated four- 
quark shape-functions v s>a (u) should not lead to the conclusion that these contributions 
are suppressed. In fact, the same factor of g 2 is also present in the tetra- local piece of 
the two- derivative matrix element (|32|). which is g 2 h v A+A+h v in light-cone gauge, and 
here the g 2 is conventionally assumed to be normalized at the non-perturbative strong 
interaction scale. 

The B decay distributions are most conveniently expressed in terms of the scalar 
components of the hadronic tensor. In the following we give the result for the currents 
relevant to semi-leptonic B — > X u iv decay and the radiative B — > Xgj decay in the 
convention specified by (j3J). For the semi-leptonic decay 

Ti = 7 ^(l- 75 ), 1^ = 7^(1-75), (62) 



and we find 



W 1 



n+p 



duo 5(n^p — uj) 



1 + 



n + p I 



Skin(^) + C mSLg (m b /jj)s mSLg (uo) 



2m b 



uS(u)-t(u) u 3 (u) + u a (u) 2(v s (u) - v a (u)) 



Wo 



Wa 



(n+p)' 



m b n+p 

— f duo 5(n^p — lu) S(uo) 

n + p J 

du 8(n-p — uj) t(u) , 
ujS{uj) 



n+p 



(n+p)' 



du 5(n_p — uj) 



t(uj) t(uj) + u a (uj) - 2v a (uj) - 
n+p 

(n + p + n_p)/2 and p 2 



(63) 



The result is given in the frame p+ = 0, where vp = [n+p + n_p)/'Z and p" = n + pn_p. 
This can be used to convert the expressions to hadronic variables, see (@J). We recall 
that the expressions are valid in the kinematic region where n + p ~ m& and n_p ~ Aqcd- 
Note that only W\ has a leading-power term. 

Turning to the radiative decay B — > X s 7, we note that the photon momentum is 
q = E~ i n+ in the frame p± = —q± = 0. Neglecting terms proportional to g M which vanish 
when the hadronic tensor is contracted with the photon polarization vector, we can take 
the current to be 



r T 



^+,7^(1-75), 



j[7l^ + ](l + 75). 



(64) 



It is simpler here not to decompose the hadronic tensor into scalar functions. Evaluating 
the traces in (jH^j) we obtain 

s kin (u) + C mag (m b /jj,)s mag (uj) 



f + ef ) / duo 5(n_p - uo) { S(uo) + 



2m h 
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uS(u) - t(u) _ Ugjuj) + Ugjuj) _ 2{v s {uj) - Vgjuj)) 
m b n + p n + p 

When the polarization of the photon is not observed, the e^f term does not contribute 
and the photon energy spectrum reads 

1 dT _ f 2n-p \ , g kin (ra_p) + C ma g(m b //x)s ma g(n_p) 

2T dE 1 ~ V m b ) [n ~ P) 2m b 

— — (t(n-p) + u s (n-p) + u a (n_p) + 2(v s {n-p) - v a (n-p)) (66) 
nib v ' 

with n^p = nib — 2E 1 . This result is valid at tree level, and in the approximation where 
the four-quark operators in the weak effective Hamiltonian are neglected. 3 

The hadronic tensor involves a power correction {s^ n {u)+C^ g {mb/ ij)s m&g {uj)) / {2m b ) 

(2) 

from the insertion of the l/m& corrections to the HQET Lagrangian, 
it is conventional to evaluate the soft matrix element with the leading-power HQET 
Lagrangian. This is advantageous in applications of HQET, where use is made of the 
heavy-quark spin-flavour symmetries. In the present case it is more convenient to not 
treat the power corrections to the HQET Lagrangian in the interaction picture. Then 
the above-mentioned term should be omitted, but the matrix elements (B\ . . . \ B) are 
evaluated with the exact HQET Lagrangian including £^q ET . In this picture the HQET 
matrix elements have a (small) nib dependence, but this is not an issue as long as 
all corrections up to a required order are included. It is further useful to regard the 
shape-functions as functions of the hadronic variable n_P = n_p + Mb — nib {P is the 
total momentum of the hadronic final state), i.e. for any function f(u) above we define 
f(u)) = f(u) + Mb — nib), such that a physical spectrum such as (jHHjl is expressed in 
terms of f(n_P). Note that Mb — nib is not the HQET parameter A, but the difference 
between the physical meson mass and the heavy quark mass. The distinction is relevant 
at the level of power corrections to the spectrum. The advantage of taking the matrix 
elements with respect to the exact HQET Lagrangian is that the support of the functions 
/ is then from to oo. 4 

Comparing our result with previous work ^21 El Ej , we find agreement for the pho- 
ton energy spectrum in the radiative decay and the hadronic invariant mass spectrum 
in the semi-leptonic decay [T3] , provided we neglect the effect from the four-quark shape- 
functions v Sja (uj). However, our general result (jB3j) and the lepton energy spectrum in the 
semi-leptonic decay following from this result differ from ^3] even when the four-quark 
shape-functions are neglected. The short-distance expansion of the hadronic tensor was 
obtained in previous work by direct matching of QCD to heavy quark effective theory 

3 Thcre is no tree level contribution from four-quark operators at leading power. At order l/mj a non- 
zero contribution arises from soft gluons attached to the charm quark loop. The degree of non-locality 
of these terms depends on whether an expansion in l/m c is performed. 

4 The upper limit is infinity as a consequence of factorization, which removes the physical upper 
limit and replaces it by a cut-off. However, dimensional regularization does not provide a dimensionful 
cut-off. 
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without the intermediate use of soft-collinear effective theory. In the tree approxima- 
tion the two approaches should give the same result. However, some of the results in 
[T2*l IT3j are obtained by an expansion in transverse momentum of a hadronic tensor that 
effectively includes the integration over neutrino momentum. This can lead to incorrect 
results, because the transverse momentum relevant to the expansion is the transverse 
momentum of partons relative to the jet, not the jet and the neutrino. The calculations 
of the photon energy spectrum in the radiative decay and of the hadronic invariant mass 
spectrum in the semi-leptonic decay are not affected by this problem. 



7 Remarks on factorization beyond tree level 

We have just seen that at tree level many simplifications take place that reduce the degree 
of non- locality of shape-functions appearing in the l/m& corrections. In this section we 
investigate whether this simplification persists beyond tree level. Since the conclusion 
will be negative, it suffices to illustrate this point for the case of abelian gauge fields. 

The manipulations below rely on the analogue of the QED Ward identity for the lead- 
ing power (abelian) SCET Lagrangian after the collinear field redefinition that decouples 
soft and collinear fields. Defining 

j + - 

J - = ^S^TSS^- (67 > 

the Ward identity for jet-functions reads 

n_0 w <O|T{j + (*)x(*i) • • • x(x m )}\0) + d^(0\T{ J^{z)x{x x ) . . . x(x m )}|0> 
+n + d {z) (0\T{j4z)x(x 1 ) . . . X (x m )}\0) 

m 

= - J2(-l) P ^\z - x k ) (0\T{ X ( Xl ) . . . x(a: m )}|0), (68) 
fc=i 

where P = if x = £ an d P = 1 if X — £■ This generalizes to insertions of composite 
operators by taking the local limit of a product of fundamental fields. (For derivative 
operators some of the delta- functions terms will then be derivatives of delta-functions.) 
Our main application of the Ward identity will be in the form 

d i zf(z + ,z ± )n„d {z) (0\T{(^O(z)x(xi)...x(x m )}\0) 

= J d 4 z[d^f(z + ,z ± )} (0\T{r ± (z) X (x 1 )... X (x m )}\0) 
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m 

- E(-!) P f( x k+, x ki _) (o|r{x(xO . . . x Wjlo). 

fc=l 



(69) 



1. Degree of multi-locality. The highest degree of non- locality in the product of soft 
fields at order l/m& comes from the double insertions of and . The double 

insertion of the mixed collinear-soft quark Lagrangian is tetra-local even at tree 
level. On the other hand 

(B\T{(h v TlW^){m c T j h v ){0) ~ij d 4 z 1 Cf\z 1 )ij d 4 z 2 cf\z 2 )}\B), (70) 

while superficially tetra-local, actually reduces to a tri-local term at tree level. Does this 
hold to all orders? 

In the abelian theory the collinear Wilson lines in Cp drop out. We use ()27|) and 
integrate by parts to write (after the collinear-soft decoupling field redefinition) 

ijafzCf{z) = J d'zz^ (Y^iD; ± Y])(z_)in„d {z) (^)(z)- (71) 
Applying the Ward identity twice to (J7UJ), we can cast this contribution into the form 

-\jd% J d^<B|(MTl)(z_)^ 

x(0\T{(Wk){x) a (£W e ) (O^J^i) r ± (z 2 )}\0) + 5-function terms, (72) 

where "5-function terms" refers to contributions from the Ward identity which are only 
bi- or tri-local. In this form it is evident that the tetra-local term vanishes at tree level, 
because the jet-function with two insertions of the transverse current integrated over the 
transverse positions vanishes in this approximation. In momentum space, the derivative 
in the definition of must be proportional to the transverse component of a collinear 
vector, but in the frame where p± = there is no such vector. By the same line of 
reasoning the jet-function cannot be expected to be zero beyond tree level. When the 
insertion of the current appears inside a collinear loop diagram, the transverse derivative 
produces factors of transverse loop momentum and an even number of such factors results 
in a term after integration. We therefore conclude that beyond tree level tetra-local 
terms appear also from the double insertions of the sub-leading collinear Lagrangian, 
and the complexity of the expression is not reduced. 

2. Transformation to the covariant derivative basis. The SCET Lagrangian has many 
terms in which the combination x^n^gF \ u (x _) occurs. These can be removed by apply- 
ing the identity (J2*7j) (see also (J7TJ)). At tree level we have seen that the various terms 
organize themselves somewhat miraculously into covariant derivatives (YHD^Y)(z) de- 
spite the fact that this object represents two separate entities (see the two terms on the 
right hand side of (J2HI))- For this reason we have chosen the covariant derivative basis 
for the soft matrix elements in Section 0] 
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Using the Ward identity it is straightforward to show that at order A, the two terms 
j[ 1)j J^ and J^J^iC^ can be combined into an expression containing a soft matrix 
element with a single covariant derivative (as in (|31|)). At order A 2 , we find (using the 
short-hand notation of Section EJ) 



jftj<o) + j(Dt j( o) + j( o) t j(o) ijC w + j( o )t j( ) i£ m . C (i) 



1 

4 tr 



p r f ^r 



(II) + (VII) + (XI) + (XII) 

i J d 4 xe tpx x 

'd 4 Zl J d 4 z 2 (S|(^y)(x_)(yt^y)( Zl _)(yt^y)( Z2 _)(yt^)( )| J B) 

x(0\T{(ZW c )(0)(W}O(x) Jl(^)}\o) 
2i [ d'z (B\CKY)(x_)(Y^tD s ± ) 2 Y)(z_)(Ylh v )(0)\B) 



x 



(o|T{(eiv c )(o)(wje)(*) (e 



in+D r 2 



*)(*)}|0> 



(73) 



where the ellipses denote further terms involving d±A±. Here we recognize the tetra-local 
shape-function with two covariant derivatives defined in (|32j). As mentioned above and 
seen explicitly in Section H this term is multiplied by a jet-function that is zero at tree 
level in the frame p± = 0. The second term does not vanish at tree level and reproduces 
the corresponding terms in Section 

It therefore appears that in general the structure of terms that follows from the cur- 
rents and Lagrangian insertions of SCET can be systematically simplified by eliminating 
all the terms that involve factors of position from the multipole expansion. Application 
of the Ward identities then organizes the terms into covariant derivatives as displayed 
above. While we did not prove this statement in general, and in particular for the 
non-abelian case, we may note that when the hybrid momentum-position space formu- 
lation of SCET is employed 6\, we are led directly to an expansion in terms of covariant 
derivatives. 



3. Convolutions with "effective" shape- functions. In the tree approximation the had- 
ronic tensor (foT?|) can be written in terms of the "effective" shape-functions defined in 
(I60|) that depend only on a single variable. This seems to be a simplification, since the 
information about hadronic matrix elements is encoded in a set of single- rather than 
multi-variable functions. 

This can always be done for contributions to the hadronic tensor from two-body cur- 
rents in SCET. In these cases n + p is the only large collinear momentum on which the 
jet-function can depend, and Lorentz invariance and dimensional analysis imply that the 
dependence on n + p can be factored out except for powers of logarithms of n + p/fi. The 
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dependence on the remaining external momentum n_p can always be formally decou- 
pled by introducing / du 5(n_p — to) to define effective shape-functions that depend on 
only one variable to. The so-defined effective shape-functions are in fact convolutions of 
jet-functions and soft matrix elements as exhibited in ()60|) in the tree approximation. 
However, beyond this approximation there are contributions from three- and four-body 
currents as given in Section El The jet-functions are then convolutions of a number of 
large collinear momentum components of the effective vertices and the dependence on 
n + p can no longer be factored, but arises after the convolution with hard coefficient 
functions. Even restricting attention to the contributions from two-body currents, the 
definitions of effective shape-functions such u a>a (u) and v s>a (uj) in (|6T)|) are useful only 
in the tree approximation, because the weight functions (jet-functions) J2,3 change at 
order a s . In any case, since the number of functions exceeds the number of observables, 
one must return to their definitions in terms of multi-local hadronic matrix elements in 
order to obtain an estimate from a particular hadronic model. 

8 Estimate of four-quark contributions 

In this section we shall make a few remarks concerning the phenomenological implications 
of our results. We do not perform an exhaustive analysis, rather we restrict ourselves to 
the new contribution from the four-quark shape-functions, which has not been studied 
in [T2J |^j . The contribution we are going to consider involves B meson matrix elements 
(B\h v (x^)h v (0)q(zi_)q(z 2 -)\B) of two heavy quark and two soft light quark fields, all 
at different positions on the light-cone. An interesting aspect of this contribution is 
that the soft quark flavour originates from the weak decay vertex. For semi-leptonic 
B — > X u lv decay it has to be a u quark. The matrix element is therefore expected to be 
different for the decay of a charged B meson (where the flavour of the spectator quark 
matches u) and a neutral B meson. As we have seen, this effect is of order 1/rrib in the 
decay spectra, while such differences in the total semi-leptonic rates between charged 
and neutral B mesons are at least of order 1 / m\ |T4*1 12T] . 

The general expression for the four-quark shape-function has been given in (J33|) . It 
involves 16 scalar functions, but only two of these, Vio and V13, appear in the hadronic 
tensors of the inclusive decays B — > A s 7 and B — > X u iv at tree level. These functions 
are tetra-local and their Fourier transforms depend on three variables. However, the 
result can be written in terms of two effective shape-functions v s>a (uj), see ffHU|) . 

In order to obtain a quantitative estimate of these contributions we have to employ a 
model, which at present can be only quite crude. Matrix elements of four-quark operators 
are often approximated in the "vacuum saturation approximation," where (leaving out 
spin and colour indices) 

{B\h v {xJ)u{z 2 J)u{z 1 J)h v {Q)\B) -> (B\h v (x-)u(z 2 ^)\0)(0\u(z^)h v (0)\B). (74) 

This gives Vg_i6 = 0, since for these functions the heavy and light quark are in a colour- 
octet configuration. We conclude that the tree level contributions are suppressed, and 
that we need a better model for a quantitative estimate. 



28 



To account for the suppression of the colour-octet matrix elements, we introduce a 
parameter e and write (now including colour) 



(B\h v {x-)T A u(z 2 -)u(zt-)T A h v (0)\B) 

-> e(B\h v (x^)u(z 2 -)\0)(0\u(z^)h v (0)\B). (75) 

For local four-quark operators the agreement of calculations of the lifetime difference 
of the B~ and B° mesons |22J EHj with observations suggests that \e\ cannot be larger 
than 0.1 without invoking large cancellations between different terms. Assuming ()75j1 
we can express the four-quark shape-functions in terms of the square of the B meson 
light-cone distribution amplitude. The resulting parameterizations of v Sta (ui) exhibit 
a delta- function singularity at w = — Mb, which can be traced to the fact that 
in the vacuum saturation approximation no dynamical (soft) momentum is transferred 
between the h v u and the uh v configuration. This is clearly not the case in reality, where 
soft gluons must be exchanged simply to rearrange colour. This motivates our "modified 
vacuum saturation approximation," where we allow the intermediate state to carry soft 
momentum k, and write 

(B\h v (x^T A u(z 2 ^u( Zl ^T A h v (0)\B) 

'dk + f(k + )(B\h v (x^u(z 2 _)\k}(k\u( Zl .)h v (0)\B}, (76) 



where \k) represents a soft state of momentum k and k + = n-k. f(k+) (whose inte- 
gral is normalized to 1) represents a (non-perturbative) distribution function of the k + 
momentum component of the soft state. We then write 

(k\u(yJ) a h v (xJ)p\B) = e -i(MB-m b -k + ), + _ x _) q ^( ) /3 | j B) 

« e-W*-**-^** (0\u{y--X-) a h v {0)j3\B) 



-i{M E 




(t) + 4- - term (77) 



with t = (y — x) + . The transition to the second line involves the assumption that the 
t-dependence of the matrix element is not modified by the "soft" state. This may be 
difficult to justify, but it allows us to obtain the desired matrix element in terms of the 
light-cone distribution amplitude of the B meson [2U 123 about which a few things are 
known. A second possible structure proportional to ?/_ is not given above, because it 
drops out in the following calculations. The momentum space distribution amplitude is 
defined by 

<p B+ (t) = / due-** (78) 



o 



It is now straightforward to insert (|77jl into (|75jl and to match the Dirac structure 
to those defining the general decomposition of the four-quark shape-function. We then 
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find that the two relevant effective shape-functions are given by 

v s (n^p) = -v a (n-p) = f B M B £ J dk + f(k+)< 5(n^P - k+) J — 0b+M 



2(j) B+ (n-P- k+) p r duj _ 

n_P — k + Jo n_P — u — k. 



na s f B M B 



\ 2 

A B 



ef(n-P) 1 



2A| 



f(n-P) 



rn—P 

x / dk+ /(n_P — k 
Jo 




(79) 



where P denotes the principal value and 1/X B — Jq°° dui (pB+ij^) 1^- We observe that the 
final result is expressed in terms of the hadronic variable n^P = P + = Eh — \Ph\- We can 
get an idea of the magnitude of the power correction from the four-quark shape-functions 
from ()59|). which instructs us to compare 4v s (n_p)/n + p to S(n_p) = S(n_p+M B —m b ) = 
S(ri-P). The integration over n + p replaces n + p by an average value (n+p) of order M B , 
hence 

Av s (n_p) f B M B /(n_P) 



c2 
f B 

A B 



■{! + ...}, 



(n + p)S(n_P) ' X B " (n + p) S(n_P) 1 ~ ' " ' 1 ' (80) 

where the ellipses in brackets stand for the twofold integral in (|79|). If we set the bracket 
to one and assume that / and S have similar shapes, we see that the size of the correction 
is about 
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±5%, 



(81) 



compatible with a power correction. Evidently there is a large uncertainty associated 
with this estimate, since, besides the assumptions in the construction of the model, 
a s could be anything between 0.3 (corresponding to the hard-collinear scale) and 1 
(corresponding to the soft scale) and |e| = 0.1 is just a guess. We should emphasize that 
the estimate is for the semi-leptonic decay of a E>~ . For B° the correction vanishes in 
our approximation. The distinctive feature of the four-quark contributions is therefore 
a difference in the B~ and B° spectra, which could be of order 5%. 

The effect could be more important than the estimate suggests, if it leads to a signif- 
icant distortion of the spectrum in the region of n_P ~ Aqcd- For instance, if f(n_P) 
were peaked at smaller arguments than S(n_P), as would be the case in the naive vac- 
uum saturation approximation, where the "soft" intermediate state has no momentum, 
the four-quark contribution could give a significant enhancement at small values of ri-P. 
To be definite, let us assume the following models for the various functions: 



S(cj) 



<f>B+(u) 



32lu 2 

7T 2 A 3 
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2A 
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vrA 2 



30 



-0.5 
-1 




0.8 1 

P+ (GeV) 



Figure 2: Distortion of the P+ spectrum in B~ — > X u iv decay by four-quark 
contributions assuming the model (|82|> . The solid central curve is S(P+). The 
dashed curves correspond to A = 100 MeV (long dashes) and A = 500 MeV (short 
dashes). Each pair is for e = ±0.1. 



r 



-k+/\ 



(82) 



To display the effect on a decay spectrum, it is natural to consider the distribution in 
the hadronic light-cone variable P + = Eh — \Ph\, since it is proportional to S(P + ) at 
leading power. Including the four-quark shape-functions but neglecting the other power 
corrections, we obtain for E>~ — > X.lv 



1 dT 
TdK 



S(P + ) 



10v s (n_p) + 2v a (n_p) 
3rrib 



(83) 



In Figure El we show the effect of the four-quark contribution on the spectrum of B~ — > 
X u tv decay. We choose the parameters a s = 1, f B = 200 MeV, A = 500 MeV, e = ±0.1 
and plot the spectrum for A = 100 MeV and A = 500 MeV. The solid curve is the 
spectrum without the power correction, which in the present approximation is also the 
decay spectrum for the neutral B° meson decay. For the smaller value of A the effect can 
be large, but it is concentrated at P + < 0.5 GeV as expected. Clearly, this is a rather 
model-dependent statement. Since the integral of f S)a (n_p) vanishes, the effect changes 
sign at some P + . Once a sufficiently large interval is integrated over, the integrated 
effect is not expected to exceed the estimate ()81|). However, we still expect the four- 
quark contributions to provide the largest difference between the spectra of charged and 
neutral B meson decays. 



9 Conclusions 

Using the framework of soft-collinear effective theory, we have investigated the l/m& 
corrections to inclusive heavy-to-light transitions in the endpoint region, where the heavy 
quark expansion in local operators breaks down. We find that SCET factorizes short- 
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and long-distance effects also at sub-leading power. The hadronic tensor, from which all 
decay spectra are derived, is represented as convolutions of hard coefficient functions, 
perturbative jet-functions and non-perturbative shape-functions. However, the structure 
is significantly more complicated than at leading power, because the factorization formula 
involves multiple convolutions and many shape-functions. We have given the general 
form of the effective currents and shape-functions that can appear at order l/m&. 

We then computed the hadronic tensor to order 1/mj, in the tree approximation. 
In this approximation the power corrections can be parameterized by a few "effective 
shape-functions" of a single variable P + only. Compared to earlier results, we find a new 
contribution from four-quark operators. The other contributions agree with previous 
results for some differential distributions, but disagree with others. We estimated the 
numerical effect of the four-quark shape-functions. The effect can lead to a distortion 
of the P + spectrum of B~ — > X u iv relative to the one of B° — ► X u £u. In our model 
the distortion can be significant at small P + and could be of order a few percent when 
integrated over an interval of order Aqcd- (The effect vanishes when integrated over 
all P + .) Our results suggest that it should suffice to include the l/m& corrections in 
the tree approximation in the analysis of measurements. This is fortunate, since a pro- 
hibitively large number of unknown shape-functions is expected to appear once radiative 
corrections are included. 

While this paper was being written, two articles [2111211 appeared that also address l/m& 
corrections to B decay distributions in the soft-collinear effective theory framework. The 
paper by Lee and Stewart [26J is similar in scope to the present work. Since they use 
a different representation of SCET, we find it difficult to compare the results in detail, 
but we can nonetheless make several comments. The structure of factorization as given 
by our equation agrees with their equation (119) on the point that the most general 
factorization formula contains convolutions over additional tri- and tetra-local shape- 
functions whose jet-functions vanish at tree-level. However, it appears to us that their 
basis of SCET currents at relative order A 2 cannot be complete because it does not 
include four-body operators, which implies that the set of power-suppressed jet-functions 
quoted in their equation (119) is also incomplete. On the other hand, Lee and Stewart 
observe that q s in can be replaced by r/ + r/_g s /4, which implies a simplification of the 
general decomposition of the four-quark matrix element (|3*3*|) . We find that only 6 out 
of the 16 invariant functions defined in (|3*3*)) can appear to any order in the perturbative 
expansion. (Two of the 8 functions defined in j2H] can be eliminated.) We differ from 
j2S] on the phenomenological implications of our results. On the basis of power counting 
it is estimated in j2E] that the effect of the four-quark contributions is up to 180% in 
large disagreement with our estimate of 5% based on a simple model. The discrepancy 
can be attributed to the absence of the colour suppression factor e in their estimate 
and numerical factors. Bosch, Neubert, and Paz (2ZI restrict their analysis to the tree 
approximation, but use the position space SCET formalism as we do. Our Section |U] has 
significant overlap with their paper and the tree level results are in complete agreement 
(as is most evident by comparing our (J55j) to their equation (28)). Our (effective) shape- 
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functions can be obtained from theirs by making the replacements 



S(u) -> S(-u) 
, N s kiQ (-u) + C mag (m b / fi)s mag (-u) 

t(u) -> 

-> u a (-uj) - 2v a (-uj). (84) 

These authors perform a more extensive phenomenological study of decay spectra, but 
they do not consider the four-quark contributions explicitly. Rather, they absorb v s>a into 
a redefinition of u Sja (their u, v) . While this is technically possible, it is also misleading, 
because it hides the fact that the four-quark contributions are different for charged and 
neutral B meson decay. In a recent paper [28J Neubert estimated the numerical effect of 
the four-quark shape-functions in various models, including the naive vacuum saturation 
ansatz (f(k + ) = S(k + )). When integrated over a sufficiently large region in P + , his result 
is in qualitative agreement with ours. 
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